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ABSTRACT 

We have carried out a systematic, homogeneous comparison of optical and near-infrared dis- 
persions. Our magnitude-limited sample of early-type galaxies in the Fornax cluster comprises 
1 1 elliptical and 1 1 lenticular galaxies more luminous than Mb = — 17. We were able to de- 
termine the central dispersions based on the near-infrared CO absorption band head for 19 
of those galaxies. The velocity dispersions range from less than 70 km s _1 to over 400 km 
s~ 1 . We compare our near-infrared velocity dispersions to the optical dispersions measured by 
Kuntschner (2000). Contrary to previous studies, we find a one-to-one correspondence with 
a median fractional difference of 6.4%. We examine the correlation between the relative dust 
mass and the fractional difference of the velocity dispersions, but find no significant trend. 
Our results suggest that early-type galaxies are largely optically thin, which is consistent with 
recent Herschel observations. 

Key words: galaxies: clusters: individual: Fornax - galaxies: elliptical and lenticular, cD - 
galaxies: kinematics and dynamics 



1 INTRODUCTION 

The dust content of early-type galaxies (ellipticals and lenticulars), 
and the effects on galaxy parameters inferred from optical observa- 
tions, is still unclear and a hot topic of debate. Dust in early-type 
galaxies was first observed in the form of dust lanes and patches 
(Hawarden et al. 1981; Ebneter & Balick 1985; Veron-Cetty & 
Veron 1988). Later HST observations revealed that dust extinction 
features exist in a large fraction of early-type galaxies (van Dokkum 
& Franx 1995). Nonetheless, the optical dust features could not ex- 
plain the high far-infrared (FIR) fluxes found by Goudfrooij & de 
Jong (1995) (based on Infrared Astronomical Satellite (IRAS) ob- 
servations). Their results implied dust masses exceeding the values 
from optical extinction by nearly an order of magnitude, indicating 
that the major part of the dust is diffusely distributed. Furthermore, 
the dust mass estimates based on IRAS flux densities are only a 
lower limit for the true dust masses, because IRAS is not sensi- 
tive to cold dust. FIR observations of early-type galaxies selected 
from the Infrared Space Observatory (ISO) archive found that the 
colder dust component dominates the total dust mass, which is typ- 
ically more than 10 times larger than the dust masses previously 
estimated using IRAS observations (Temi et al. 2004). Recent FIR 
Spitzer observations of elliptical galaxies show evidence of diffuse 
dust (Temi et al. 2007) and ground-based submillimetre continuum 
observations with SCUBA have already revealed that galaxies in- 

* Based on new observations collected at the European Southern Observa- 
tory, Cerro Paranal, Chile; ESO program 082.B-0897A. 



deed contain large amounts of cold dust, but that the submm emis- 
sion of some of the elliptical galaxies may be synchrotron rather 
than dust emission (Vlahakis et al. 2005). Nevertheless, the Her- 
schel Space Observatory is crucial for detecting the missing cold 
dust component and, as a consequence, for making accurate dust 
estimates in ellipticals (Boselli et al. 2010). Science Demonstration 
Phase results for the Herschel Virgo Cluster Survey (Davies et al. 
2010) did not confirm the diffuse dust thesis in elliptical galaxies: 
Clemens et al. (2010) found no detection for passively evolving 
early-type galaxies and Baes et al. (2010) did not find evidence for 
a diffuse dust component in M87, explaining the FIR emission by 
a synchrotron model. 

The influence of dust on the optical and NIR photometry of 
early-type galaxies was studied by Michard (2005) and models 
to predict the effect of dust were developed by Witt et al. (1992) 
and Wise & Silva (1996, 1997). Baes & Dejonghe (2000, 2001, 
2002) and Baes et al. (2000) were the first to include the influ- 
ence of dust on the observed stellar kinematics in their models. 
They showed that dust may bias optical observations through ab- 
sorption and scattering, thus influencing the photometric and kine- 
matic data. Such effects would have widespread ramifications for 
studies of early-type galaxies, as velocity dispersions play a sig- 
nificant role in tracking the mass evolution of early-type galaxies 
(van der Marel & van Dokkum 2007), and appear in empirical re- 
lationships such as the the Mbh-o" relation (Ferrarese & Merritt 
2000; Gebhardt et al. 2000; Graham 2008; Gultekin et al. 2009; 
Kormendy & Bender 2009; Graham et al. 2010), the Faber- Jackson 
relation and the Fundamental Plane relation (Djorgovski & Davis 
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1987; Dressier et al. 1987; Guzman et al. 1993; Pahre et al. 1995; 
Prugniel & Simien 1996; Graham & Colless 1997; Gavazzi et al. 
1999; Bernardi et al. 2003; Michard & Prugniel 2004; Desroches 
et al. 2007; Fraix-Burnet et al. 2010; Ribeiro & Dantas 2010). 

One way to investigate the effect of dust on the kinematics is 
to compare the optical and NIR velocity dispersions. The K-band 
is the optimal region to perform a study on NIR velocity disper- 
sions: the extinction by dust in the K-band is only 7% of that in 
the B-band (where dust is very opaque) and the wavelengths are 
short enough to evade dilution of the stellar continuum by hot dust 
(Gaffney et al. 1995). In the K-band, one can use the 2.29 /im 
(2—0) 12 CO absorption band head as a kinematical tracer. It is the 
strongest absorption feature in galactic spectra in the 1—3 fim range 
and increases in strength with decreasing effective temperature or 
increasing radius of the underlying stars (Silge & Gebhardt 2003). 
This late-type star feature can be used to probe the kinematics of 
the red stellar population of galaxies (Gaffney et al. 1993) and is 
intrinsically sharp and deep, spectrally isolated from other strong 
absorption and emission features and is located in a dark part of the 
infrared sky spectrum (Lester & Gaffney 1994). The advent of effi- 
cient NIR spectrographs and detectors on large telescopes has made 
the CO absorption band a generally available tool for kinematic 
studies (e.g. Silge et al. 2005; Nowak et al. 2007, 2008; Lyubenova 
et al. 2008). 

Comparative studies about NIR and optical velocity disper- 
sions were performed in the past. Silge & Gebhardt (2003) were 
the first to perform a systematic study, but their sample was biased 
to lenticulars (it contained 25 galaxies, of which 7 ellipticals and 18 
lenticulars) and was based on inhomogeneous optical velocity dis- 
persions, obtained from different papers with different instruments, 
methods and extraction windows. They found that the velocity dis- 
persion decreases with wavelength, opposite to the theoretical ex- 
pectations of Baes & Dejonghe (2002). Silva et al. (2008) studied 
the stellar populations in early-type galaxies, using the strong spec- 
tral features near 2.2 /im. They found a one-to-one correspondence 
between the optical and NIR velocity dispersions, based on a ho- 
mogeneous set of 4 lenticular and 7 elliptical galaxies. 10 out of the 
1 1 galaxies are part of our new, extended sample, in which we ex- 
cluded NGC 1344, because this galaxy was not included in the sam- 
ple of Kuntschner (2000). Romberg & Fischer (2010) inspected the 
NIR and optical velocity dispersions in infrared-luminous mergers. 
Their reference sample of elliptical galaxies did not show much ev- 
idence for a CT-discrepancy, but they found a strong correlation with 
IR-luminosity and dust masses in merging galaxies. 

To investigate the disagreement between the different veloc- 
ity dispersion studies, we embarked on a project, for the first time 
using a complete and well-balanced sample of 22 early-type galax- 
ies (comprising 11 ellipticals and 11 lenticulars). All galaxies are 
members of the Fornax cluster, the nearest galaxy cluster after the 
Virgo cluster. Fornax is considerably more compact and regular in 
shape than Virgo, doubling the central density of galaxies but hav- 
ing a total mass of nearly an order of magnitude lower. This makes 
the Fornax cluster a good representative of the groups and poor 
clusters in which most galaxies in the universe reside (Jordan et al. 
2007). Nevertheless, the main reason for the choice of our sample 
was the availability of uniform spectroscopic optical data, obtained 
by Kuntschner (2000). In his paper, velocity dispersion, age, metal - 
licity and line strength are discussed for all galaxies in the sample. 

This paper is organized as follows: § 2 presents the data and 
the data reduction, § 3 presents the comparison of the optical and 
IR velocity dispersions and the relation with the dust mass. In § 4 
our conclusions are presented. 



2 OBSERVATIONS AND DATA REDUCTION 

2.1 Data sample 

This study is based on the complete magnitude-limited sample of 
Kuntschner (2000). This sample, presented in Table 1, has been se- 
lected from the catalogue of Fornax galaxies of Ferguson (1989), 
in order to obtain a complete sample down to B T = 14.2 or 
Mb = —17 and contains 11 elliptical and 11 lenticular galaxies. 
Recent measurements of Bt were taken from the NASA/IPAC Ex- 
tragalactic Database. 

2.2 Observations 

The observations were performed between 20 October 2008 and 26 
January 2009 with the VLT using the SW arm of ISAAC (Moor- 
wood et al. 1998) in spectroscopic medium resolution mode (SWS- 
MR). The characteristics of the detector and the instrumental set-up 
are given in Table 2. 

We used the nod-on-slit mode (double subtraction technique). 
Every observational sequence was started with the galaxy centered 
on the slit near one end and an individual spectrum was taken. The 
galaxy was then moved 60" or 90" towards the other end of the 
slit and two more integrations were executed. Then the galaxy was 
placed again at the original slit position where another spectrum 
was taken. This ABBA sequence was repeated a number of times, 
resulting in multiple individual two-dimensional spectroscopic im- 
ages. In order to remove properly the sky background for the obser- 
vations of big and bright galaxies (i.e. for NGC 1399, NGC 1316, 
NGC 1380, NGC 1404 and NGC 1427 in particular), we observed 
separate sky spectra. That way, only half of the frames contained 
the galaxy spectrum, resulting in an observing sequence of the form 
OSSO (0=object integration, S=sky integration). Total on source 
exposure times were between 900 and 2400s [see Table 1 for a de- 
tailed listing, with '(1/2)' denoting OSSO observations]. The seeing 
was generally sub-arcsecond, reaching as low as ~ 0.5" in the best 
cases. 

To be able to remove the telluric lines, we observed B-type 
standard stars with the same observing sequence we have used for 
our science targets. These hot stars do not have any spectral features 
in the wavelength interval we are considering. Moreover, their con- 
tinuum in the if-band is well approximated by the Rayleigh-Jeans 
part of the blackbody spectrum (Silva et al. 2008). To remove the 
stellar signature, we divide the spectrum by a template of the corre- 
sponding type of star (Pickles 1998). All the remaining variation in 
the spectra is due to observational features and is caused by the tel- 
luric lines and the instrument response, so, by dividing the galaxy 
spectra by this residual template spectrum, we remove the telluric 
lines and the detector signature. Because of the variability of the 
NIR sky, we impose the conditions that the difference in observing 
time between the target and the standard star has to be less then 2 
hours and that the difference in airmass between the observations 
has to be less than 0.2. 

2.3 Data reduction 

The basic data reduction steps were performed with the ISAAC 
pipeline and MIDAS. The sky subtraction, bias subtraction, flat 
fielding and wavelength calibration (based on the OH-lines) were 
done by the ISAAC pipeline. The removal of telluric lines (using 
our standard star observations) was done with MIDAS. If possi- 
ble and appropriate, we increased the S/N of the telluric profile by 
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Table 1. Galaxy sample. 
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Notes. -Morphological type and optical velocity dispersions taken from Kuntschner (2000), B T from NASA/IPAC Extragalactic Database. Empirical S/N 

estimates were determined for the extracted spectra as described in §2.3. 

approximately doubling the original step size. To illustrate the ef- 
fect of the different data reduction steps, we refer to Fig. 2, where 
a central 3.7" extraction of NCG 1381 is shown for the different 
data reduction steps: the ISAAC pipeline result, after removing the 
telluric lines and after rebinning to the 2.43A/pix wavelength in- 
crement. For each extracted spectrum, we have derived an empiri- 
cal S/N following the method described by Stoehr et al. (2007). The 
resulting S/N per rebinned element ranges from ~ 7 for the faintest 
galaxies to more than 100 for the brightest ones and is listed in Ta- 
ble 1. 

The instrumental set-up, as described in Table 2, gives a 
FWHM resolution of approximately 107 km s~ , producing an in- 
strumental contribution to the dispersion of ~ 45 km s _1 . 

2.4 Velocity dispersion determination 

The sharp blue edge of the CO band head allows us to mea- 
sure the kinematics accurately (Silge & Gebhardt 2003). There are 
several techniques to obtain the internal kinematical information, 
e.g. the Fourier correlation quotient (FCQ) method, developed by 
Bender et al. (1994) and used by Kuntschner (2000). We use the 
pPXF technique developed by Cappellari & Emsellem (2004). This 
method rebins the spectrum logarithmically and fits it directly in 
pixel space. We correct the template continuum shape using ad- 
ditive second degree Legendre polynomials and use a pure Gaus- 
sian to model the LOSVD. An alternative LOSVD is inspected in 
§ 3.2, where Gauss-Hermite coefficients fe 3 and h$ are included 
to model the LOSVD. We apply pPXF to a wavelength range of 
~ 2.23 — 2.345 /im, fitting not only the CO band head, but also 
including a part of the spectrum bluewards of it. To illustrate the 
effect of the velocity dispersion on the CO band head, we refer to 
Fig. 2 of Silge & Gebhardt (2003). 

The choice of the template stars, which we obtained from the 



Table 2. The instrumental set-up. 



Telescope 

Dates 

Instrument 



VLT UTl/Antu 

20 October 2008 - 26 January 2009 
ISAAC 



Spectral range 
Grating 
Dispersion 
Resolution (FWHM) 
Spatial Scale 
Slit Width 
Detector 
Gain 
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Seeing 
Pixel Size 
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1 arcsec 
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4.6 e" ADU" 1 
~ 10 e" 
< 1" 
18.5/im 



averaging different observations of telluric stars. For NGC 1380, 
telluric lines have not been removed, because of a problem with the 
wavelength calibration of the standard star. For all spectra we had 
some problem around 2.3 /im to correct a telluric absorption line 
with a very sharpe edge. As a result, a spike-like emission feature 
remains in the corrected spectra. We choose not to remove cosmic 
rays, but will ignore any contaminated pixel when fitting the stel- 
lar template. An average of the 25 central rows was taken to obtain 
an effective aperture of 1 x 3.7 arcsec 2 , choosing a slit width of 1 
arcsec to obtain the needed spectral resolution and approximating 
the spatial width of the 2.3 x 3.85 arcsec optical aperture used 
by Kuntschner (2000). We do not increase the extraction width to 
obtain a better S/N, because the velocity dispersion can be strongly 
aperture-dependent (as shown in Fig. 1), and it is important to be 
well-matched to Kuntschner (2000). As a final step, the spectra 
were rebinned to a common wavelength increment (2.43 A/pix), 
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Figure 1. Central velocity dispersion of IC 2006 as a function of the extrac- 
tion width. The dashed line shows the standard extraction aperture of 3.7 
arcsecs. 
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Figure 2. Data reduction steps applied to a central 3.7" extraction of 
NGC 1381 are shown (from top to bottom) before corrections for telluric 
absorption (the pipeline result), after correction for telluric absorption and 
instrument sensitivity and after rebinning to ~2.43A/pix. The sharp emis- 
sion feature around 2.3 /im is a remnant of the telluric correction.No correc- 
tion for redshift has been applied. The middle and the bottom spectra have 
been separated by an arbitrary additive offset; otherwise, they would be on 
top of each other. 

GNIRS and NIFS libraries (Winge et al. 2009), is an important as- 
pect when determining the velocity dispersion: Silge & Gebhardt 
(2003) showed that the equivalent width of the template used for 
the fit is important, not the details of the spectral type. To account 
for the spectral resolution differences between ISAAC and the tem- 
plate libraries, we have logarithmically rebinned the latter to the 
same velocity scale as the ISAAC spectra, under the assumption 
that the shape of the instrumental spectral profiles can be well ap- 
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Figure 3. Dispersion measured by the pPXF method for IC 1963 as a 
function of the equivalent width of the input template star. The two filled 
squares are the results for an average template of K giants (EWco = 8.99, 
ao = 48.6 ± 5.6) and M giants (EW C0 = 15.01, cr = 32.9 ± 4.7). The 
asterisks represent the EWs of the galaxies in the sample. 

proximated by a Gaussian. Fig. 3 presents the measured central ve- 
locity dispersion as a function of the equivalent width of the input 
template, with equivalent widths ranging from ~ 1.5 to ~ 15. One 
concludes that the larger the EW of the CO band of the template, 
the lower the resulting velocity dispersion required to reproduce 
the galaxy spectrum. The relative values of the velocity dispersion 
are only comparable, if a fixed template is used to generate the fits. 
The filled squares are the results for an average template of K giants 
(HD206067, HD218594, HD39425 and HD4730 from the GNIRS 
library) and M giants (HD27796, HD30354, HD23574, HD234791 
from the NIFS library). To be consistent with Kuntschner (2000), 
we choose the average K giant spectrum to fit the galaxy spectra 
in a homogeneous way. An alternative template scheme is explored 
in § 3.2, where velocity dispersions have been determined using all 
the templates of Fig. 3 as input templates for pPXF. The asterisks 
in Fig. 3 represent the EWs of the galaxies in the sample, given in 
Table 1 . This shows that the stellar library spans the range of the 
observed CO strength in the galaxies. 

Table 1 presents the resulting NIR velocity dispersion unir 
and the optical velocity dispersions er opt taken from Kuntschner 
(2000). The optical velocity dispersions have been used as the ini- 
tial value in the pPXF method. After trying several approaches, we 
have used a bootstrap method to estimate the uncertainties on the 
velocity dispersions, resampling the residuals of the initial fit. For 
ESO359-G02, ES0358-G25 and ESO358-G50 we did not get a 
reliable fit, as a consequence of the low S/N for these galaxies and 
the low expected velocity dispersions (based on the optical values). 
We decided to exclude these galaxies from the sample. Afterall, 
their optical velocity dispersions are not reliable either; they suf- 
fered an instrumental broadening of ~ 105 km s _1 , which intro- 
duces optical systematic errors for a < 90 km s _1 [as stated in 
Kuntschner (2000)]. 

The spectral fits are presented in Figs. 4 and 5: the thick 
smooth line is the resulting fit, the thin noisy line represents the 
galaxy spectrum. Cosmic rays, sharp features and poorly removed 
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telluric features were flagged as bad pixels and were not used by 
the pPXF method to generate the fit. If we estimate the error on the 
fit by the dispersion of the residuals, we obtain x 2 /DOF < 1 for 
all galaxies, which indicates that we obtain good fits. 



3 RESULTS AND DISCUSSION 

3.1 Comparing NIR velocity dispersions 

The line strength study of Silva et al. (2008) had 4 SO and 6 E galax- 
ies in common with our galaxy sample [of which 8 are observed 
with ISAAC and 2 (NGC 1316 and NGC 1399) with SINFONI at 
the VLT]. For their ISAAC observations in SWS-MR mode, they 
also used a 120" x 1" slit and had total exposure times of 3000s or 
3200s. This resulted in extracted spectra with a S/N ranging from 
48 to 280, based on the S/N determination of Stoehr et al. (2007). 
However, they applied a procedure derived from the processing of 
optical long-slit spectra and did not use a simple AB subtraction, 
because this could cause inaccurate dark and background correc- 
tion in the lower surface brightness parts of the luminosity profiles. 
They periodically observed velocity template standard stars, cover- 
ing the range K5 III to Ml III. Those templates were fitted to the 
galaxy spectra applying the pPXF method, extracting the NIR ve- 
locity dispersions ONiR.siiva- Unfortunately, no uncertainties on the 
velocity dispersions were given. They used 1" x l/8R e for the 
extraction window, which results in spatial widths ranging from 
0.15 to 0.61 arcsec (extracted R e for NGC 1404 equals 2.9 arcsecs, 
kindly provided by Dr. Kuntschner, private communication). For 
the SINFONI observations, Silva et al. (2008) extracted a slit of 
2" x 3" from the data cube. 

In Fig. 1 it is shown that the velocity dispersion can depend on 
the extraction width. In order to compare our measurements with 
Silva et al. (2008), we have redetermined the velocity dispersion 
for the 10 galaxies in common, adopting an extraction width of 
1" x l/8.Re (1" x 3" for NGC 1316 and NGC 1399) and fitting 
the average K giant template to the galaxy spectra. The measure- 
ments are denoted as <7NiR,i/8Re. The results are presented in Fig. 
6, where the derived velocity dispersions aNiR,i/8Re are compared 
to the values <TNiR,siiva, derived by Silva et al. (2008). In this and 
following figures, true elliptical galaxies are represented as open 
triangles, lenticular galaxies as filled circles. However, to distin- 
guish the 2 galaxies observed with SINFONI by Silva et al. (2008), 
we have marked those galaxies with asterisks in Fig. 6. We have 
used the uncertainties on onir as an estimate for the uncertainties 
on <7nir,i/8 Re - The best-fitting line, which is given by 

^NiRj/sRe = (1.04 ± 0.05)a NrR , S uva - (13.78 ± 9.52), (1) 

has a slope consistent with 1 , an intercept marginally inconsistent 
with and a reduced x 2 of 1.89. Using the velocity dispersions 
obtained with the fixed extraction width of 3.7" (ctntr) instead of 
onir, i/8 Re does not change the fit significantly. Thus, our NIR results 
generally compare well to those of Silva et al. (2008) for the same 
galaxies. 

3.2 Comparing NIR and optical dispersions 

Fig. 7 presents the correlation between the optical and NIR central 
velocity dispersions for 19 galaxies of our sample. The solid line 

cr opt = (0.99 ± 0.06)ct N ir - (3.76 ± 8.72) (2) 




Figure 6. Correlation between the velocity dispersions determined with an 
extraction window of 1/8 R a and the NIR velocity dispersions from the 
literature (Silva et al. 2008). The dashed line has a slope of unity, the solid 
line is the best fit, given by equation 1 . Here and in the following plots, the 
filled circles are SO galaxies, while the open triangles are E galaxies. The 
two galaxies observed with SINFONI (Silva et al. 2008) are represented 
with an asterisk. 

shows the best fit, with a slope consistent with 1, an intercept 
consistent with and a reduced \ 2 of 1.09. The absence of a a- 
discrepancy confirms the findings by Silva et al. (2008) and fur- 
ther generalizes their results due to the statistical completeness of 
this study, but does not agree with Silge & Gebhardt (2003). Their 
best-fitting line had a slope of 1.189 ± 0.084 and an intercept of 
-8.6 ± 12.4. 

Kuntschner (2000) had difficulties in determining the velocity 
dispersions for galaxies with a opt < 70 km s _1 and noted that the 
results of ~ 50 — 60 km s _1 are only rough estimates. However, 
excluding the galaxies with <j opt < 70 km s -1 does not change the 
results significantly. 

Fig. 8 presents a histogram of the fractional difference be- 
tween NIR and optical measurements of dispersion. The median 
fractional difference between the optical and the NIR velocity dis- 
persions is 6.4%, the mean fractional difference is 3.9%. Silge & 
Gebhardt (2003) found a median difference of —11%, opposite to 
theoretical expectations (Baes & Dejonghe 2002) and to our results. 

It is not immediately clear why Silge & Gebhardt (2003) 
found that NIR velocity dispersion measurements are lower than 
optical dispersion measurements. One possible reason could be the 
different spatial range used by Silge & Gebhardt (2003): Fig. 1 
indicates that velocity dispersions can decrease for higher spatial 
ranges. The average spatial width in Silge & Gebhardt (2003) is 
~12", which may be the reason for their difference between NIR 
and optical velocity dispersions. In order to investigate this possi- 
bility, we recalculated the NIR velocity dispersions for all galaxies 
in our sample, now using an extraction window of 1" x 12". The 
best-fitting line, with a reduced \ 2 of 1.35, is 

a opt = (1.03 ± 0.07)a NIRilxl2arcscc 2 - (3.35 ± 9.16), (3) 

fully consistent with Eq. 2. Fig. 1 suggests an increase of the slope 
of ~ 4%, which is indeed what we find. Nonetheless, a ~ 20% 
effect is needed to explain the results found by Silge & Gebhardt 
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Figure 4. Rest-frame spectra of galaxies (noisy thin curves) and spectra of the average Kill stellar spectrum convolved with the derived velocity distribution 
(smooth thick curves). Cosmic rays and poorly removed telluric features were flagged as bad pixels and not used to generate the fit. No removal of telluric 
lines was done for NGC 1380. 



(2003), so we exclude that the extraction width is responsible for 
the discrepancy. 

Another possible explanation for the discrepancy could be the 
choice of the templates: Fig. 3 showed that the measured veloc- 
ity dispersion depends strongly on the equivalent width of the CO 
band of the template. Silge & Gebhardt (2003) used templates with 
equivalent widths ranging from less than 5 to over 20 A. So far, 
we have only used the average K giant template with an equiv- 
alent width of 8.99 A. In order to investigate this possibility, we 
have provided the 23 templates we have used in Fig. 3 as input 
templates for pPXF. The resulting velocity dispersions are given in 
Table 1 . For every galaxy separately, pPXF gives weights to the dif- 
ferent templates to obtain the best fit, delivering stellar population 
information. pPXF favoured HD1 13538 (KV), HD2490 (Mill) and 
HD63425B (Kill), which indicates that the light is dominated by 
cool giants and dwarfs. The best linear fit is given by 

a ovt = (0.99 ± 0.07)cr N IR,23 templates + (2.92 ± 9.05) (4) 

and has a reduced \ 2 of 1.70. It is fully consistent with equation 



2, indicating that template mismatch does not affect our velocity 
dispersion determination. 

Silge & Gebhardt (2003) used a nonparametric LOSVD. To 
investigate whether this could account for the discrepancy, we in- 
cluded the hs and hi Gauss-Hermite coefficients, using the average 
K giant template as input spectrum and the default pPXF parameter 
A to penalize both parameters. The resulting values for hs (and hi) 
are between —0.05 and 0.1 1 (—0.01 and 0.08, respectively) with a 
median value of 0.012 (0.020, respectively). The best-fitting line is 
given by 

(Top, = (0.99 ± 0.05)ct N ir,g-h - (2.40 ± 7.58) (5) 

with a reduced \ of 0.93. Again, the slope is consistent with 
1 and the intercept is consistent with 0, so we exclude that the 
parametrization of the LOSVD is accountable for the discrepancy. 

One remaining possibility for the difference between our re- 
sults and Silge & Gebhardt (2003) is that the sample selection is im- 
portant. Our Fornax sample represents only a cluster environment, 
while the Silge & Gebhardt (2003) sample ranges from cluster to 
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Figure 7. Correlation between the dispersion measured from the CO band 
head and the optical dispersion (Kuntschner 2000). The dashed line shows 
where two measurements are equal, the solid line is the best-fitting line, 
given by equation 2. 



Figure 8. Histogram showing the number of galaxies in each bin of frac- 
tional difference between infrared and optical measurements. The median 
fractional difference is 6.4%, the mean fractional difference is 3.9%. 
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isolated field galaxies and there are known correlations between en- 
vironment and velocity dispersion (e.g. Zhu et al. 2010; La Barbera 
et al. 2010). 

3.3 Correlation with IRAS dust masses 

According to Silva et al. (2008), at least two galaxies (NGC 1316 
and NGC 1380) from our sample show clear central dust features. 
We need to know how much dust there is in the galaxies to study the 
effects of dust on the observed kinematics. In this section, we derive 
the dust masses (M d ) for our galaxies based on IRAS flux densities 
at 60 and 100 fim. We apply the technique of Goudfrooij & de Jong 
(1995), using a distance D — 19.3 Mpc for all galaxies (Jordan 
et al. 2007). For 5 galaxies of our sample, both Set) and S100 were 
given. In that case both Td and Md were determined. If only S100 
was given, which was the case for 5 galaxies of our sample, we used 
the average dust temperature Td = 25. 9K of our galaxy sample as 
a representative value to compute the dust mass. Table 3 presents 
the resulting dust masses M d , T d and the IRAS flux densities for 
10 galaxies in our sample. For six galaxies (NGC 1374, NGC 1375, 
NGC 1381, NGC 1419, NGC 1427, ESO 358-G06) the IRAS Faint 
Source Catalog lists only upper limits, the three remaining galaxies 
(NGC 1373, NGC 1380A, IC 1963) are not listed in the catalog. 

Fig. 9 presents the fractional difference of the infrared and the 
optical dispersions as a function of the relative amount of dust in a 
galaxy, which is estimated in Silge & Gebhardt (2003) by the ratio 
of the IRAS dust mass to the B-band luminosity. The estimates for 
the relative amount of dust are given in Table 3. Note that the IRAS 
dust mass estimates are a lower limit for the true dust masses, be- 
cause IRAS is not sensitive to cold dust (which emits the bulk of its 
radiation longwards of 100/im.). A negligible trend is visible (con- 
firmed by the Spearman rank-order correlation coefficient equal to 
0.21). The best-fitting solid line is given by equation 

aNIR ~ a ° vl = (0.008 ± 0.056) log + (0.068 ± 0.266), 

Copt \Lb ) 

(6) 

with both the slope and the intercept consistent with and a reduced 
X 2 of 1.76. This implies that warm dust does not affect optical dis- 
persions. We cannot yet make the same conclusion for colder dust, 
but Fig. 7 indicates that the typical effect is very weak. 



4 CONCLUSIONS 

In this study, we investigate a complete magnitude-limited and un- 
biased sample of 22 early-type galaxies in the Fornax cluster and 
are able to determine the kinematics based on the 2.29/im 12 CO(2- 
0) feature for 19 of those galaxies. We related the NIR velocity 
dispersions with the optical dispersions of Kuntschner (2000) and 
found no evidence for a er-discrepancy for the ellipticals nor for the 
lenticulars. Our results agree with a previous smaller study of Silva 
et al. (2008), but not with Silge & Gebhardt (2003). We investigated 
this disagreement by providing a variety of input templates to pPXF 
with a large range of EWs, by changing the spatial width of the ex- 
traction window and by introducing Gauss-Hermite coefficients in 
the LOSVD, but we were not able to clarify this discrepancy. 

We have computed the dust masses based on IRAS flux den- 
sities for 10 galaxies of our sample and investigated the influence 
of diffuse dust on the observed kinematics, which turned out to be 
negligible. 

The one-to-one correspondence between the optical and the 
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Figure 9. Fractional difference between the infrared and optical dispersions 
as a function of the ratio of dust mass to B-band luminosity. The solid line is 
the best fit, given by Eq. 6. The Spearman rank-order correlation coefficient 
is 0.21 at a significance level of 0.56. 

NIR velocity dispersions found for this homogeneous set of early- 
type galaxies implies that velocity dispersions measured at optical 
wavelengths are reliable kinematic parameters for early-type galax- 
ies and hence that no bias is introduced in statistical relations that 
build on such dispersions (such as the Mbh-ct relation or the Funda- 
mental Plane). Combined with the simulations by Baes & Dejonghe 
(2000, 2002), it also supports the traditional point of view on the 
dust content of early-type galaxies, namely that they are virtually 
optically thin. While some observational studies hinted towards 
the existence of a substantial diffuse dust component in early-type 
galaxies (e.g. Temi et al. 2004, 2007; Leeuw et al. 2004; Vlahakis 
et al. 2005), the most recent results from the recently launched Her- 
schel Space Observatory indicate a dearth of diffuse dust in the few 
elliptical galaxies studied so far (Clemens et al. 2010; Baes et al. 
2010; Gomez et al. 2010). 

Whereas our results support this scenario, we must be careful 
for two caveats. On the one hand, our results only set limits on the 
presence of a smooth, diffusely distributed dust component, which 
one would expect if the dust has an internal origin. If the dust has an 
external origin, it is not necessarily coincident with the stellar body. 
An example is the nearby Virgo Cluster elliptical M86, where the 
dust is clearly related to stripping from a nearby spiral galaxy and is 
concentrated some 10 kpc to the south-east of the nucleus (Gomez 
et al. 2010). 

On the other hand, the comparison of optical and NIR veloc- 
ity dispersions might not be the most sensitive way to measure the 
optical thickness of early-type galaxies. The simulations of Baes 
& Dejonghe (2000, 2002) indicate an effect of a few percent only 
for optical depths of order unity. Combined with the measurement 
errors and other possible effects such as different stellar popula- 
tions dominating the kinematics at optical and NIR wavelengths, 
our results should be considered only in a statistical sense and one 
should take care not to directly interpret results on individual galax- 
ies in terms of optical thickness. A clear example is the case of 
NGC 1380: in spite of a clear dust lane in optical images (Jordan 
et al. 2007) and a significant IRAS dust mass of 7.6 x 10 5 M , the 
NIR velocity dispersion is lower than the optical dispersion. 
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Table 3. IRAS flux densities and dust characteristics. 



Galaxy 


S(60/im) 


S(100/im) 


log(M d ) 


T d 


log(M d /L B ) 




[Jy] 


[Jy] 


[Mq] 


[K] 


[M Q /L ] 


NGC 1316 


3.070 ± 0.030 


8.110 ± 1.900 


6.11 ±0.25 


27.2 ± 2.2 


-4.64 


NGC 1336 




0.260 ± 0.095 


4.73 




-4.54 


NGC 1339 


0.2298 ±0.14 


0.670 ± 0.058 


5.09 ± 0.25 


26.5 ±2.1 


-4.42 


NGC 1351 


0.090 ±0.021 


0.510 ±0.042 


5.41 ± 0.24 


22.4 ± 1.5 


-4.10 


NGC 1379 




0.140 ±0.047 


4.46 




-5.33 


NGC 1380 


1.040 ±0.042 


3.440 ±0.107 


5.88 ± 0.24 


25.6 ± 1.9 


-4.27 


NGC 1399 




0.300 ± 0.082 


4.79 




-5.48 


NGC 1404 




0.270 ± 0.056 


4.75 




-5.36 


IC2006 


0.120 ±0.015 


0.320 ± 0.047 


4.71 ±0.25 


27.2 ±2.2 


-4.92 


ES0 358-G59 




0.380 ± 0.060 


4.90 




-4.01 



A clearer picture on the optical thickness of early-type galax- 
ies will hopefully emerge in the near future when substantial num- 
bers of nearby early-type galaxies will be observed as part of sev- 
eral Herschel Key Programs, including the Herschel Virgo Clus- 
ter Survey (Davies et al. 2010) and the Herschel Reference Survey 
(Boselli et al. 2010). 
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